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Abstract 

A simple algorithm, which exploits the associativity of the BCH formula, and that 
can be generalized by iteration, extends the remarkable simplification of the Baker- 
Campbell-Hausdorff (BCH) formula, recently derived by Van-Brunt and Visser. We 
show that if [V, Y] = uX -|- vY + cl, [Y, Z] = wY + zZ + dl, and, consistently with 
the Jacobi identity, [V, Z] = mX + nY + pZ + el, then 

exp(A) exp(y) exp(Z) = exp(aA + hY + cZ + dl) 

where a, b, c and d are solutions of four equations. In particular, the Van-Brunt 
and Visser formula 

exp(A) exp(Z) = exp(aV + bZ + c[X, Z] + dl) 

extends to cases when [X, Z] contains also elements different from X and Z. Such 
a closed form of the BCH formula may have interesting applications both in math¬ 
ematics and physics. As an application, we provide the closed form of the BCH 
formula in the case of the exponentiation of the Virasoro algebra, with SL2(C) 
following as a subcase. We also determine three-dimensional subalgebras of the 
Virasoro algebra satisfying the Van-Brunt and Visser condition. It turns out that 
the exponential form of SL2(C) has a nice representation in terms of its eigenvalues 
and of the hxed points of the corresponding Mobius transformation. This may have 
applications in Uniformization theory and Conformal Field Theories. 
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1 The algorithm 


Very recently Van-Brunt and Visser [I] (see also [2] for related issues) found a remark¬ 
able relation that simplifies, in important cases, the Baker-Campbell-Hausdorff (BCH) 
formula. Namely, if X and Y are elements of a Lie algebra with commutator 

[X,Y]=uX + vY + cl , (1.1) 

with I a central element and u, v, c, complex parameters, then [T] 

exp(V)exp(V) = exp(X +V + /(M,n)[X,V]) , (1.2) 


where 


f{u,v) = f{v,u) = 


{u — v)e“+’' — (Me“ — ve^) 


(1.3) 


uv{e'^ — e'’) 

In the following we formulate an algorithm which exploits the associativity of the BCH 
formula. Set a + (3 = 1 and consider the identity 


exp(X) exp(V) exp(Z) = exp(V) exp(aV) exp(/SV) exp(Z) . 


If 


[X, Y]=uX + vY + cl , [V, Z]=wY + zZ + dI, 

then, associativity of the BCH formula and Eqs. (ll.2p(ll.4p . yield 

exp(X) exp(aV) = exp(X) , exp(/9V) expPZ) = exp(V) , 


( 1 . 4 ) 

( 1 . 5 ) 

( 1 . 6 ) 


where 


X := Qaiu, v)X + ha{u, v)Y -|- la{u, v)cl , 

Y := hp{z, w)Y + gp{z, w)Z + ^{z, w)dl , (1.7) 

with ga{u, u) := 1 -|- auf{au, v), ha{u, v) := a{l + vf{au, v)) and lain, v) := af(au, v). 
Imposing 

[X, Y]=uX + vY + dI , (1.8) 

hxes a, u, v and c. This solves the BCH problem, since, by fll.ip . fll.2p . fll.4l) and fll.61) 

exp(X)exp(V)exp(Z) = exp(X) exp(V) = exp(X + V +/(h, h)[X, V]) . (1.9) 

Note that the commutator between X and Z may contain also Y 

[X, Z] = mX + nY + pZ + el . (1-10) 

This is consistent with the Jacobi identity 

[X, [V, Z]] + [V, [Z, X]] + [Z, [X, V]] = 0 , (1.11) 

that constrains e,m,n and p by a linear system. Setting Y = XqQ and A_ := Aoa, 
A+ := Ao/3, Eq. (ll.4p includes, as a particular case, 

exp(X) exp(Z) = lim exp(X) exp(A_(5) exp(A+Q) exp(Z) . (1-12) 

Ao^-O 
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This explicitly shows that the algorithm solves the BCH problem for exp(X)exp(Z) in 
some of the cases when [X, Z] includes elements of the algebra different from X and Z. 
The complete classification of the commutator algebras, leading to the closed form fll.bp 
of the BCH formula, is investigated in [3]. The algorithm has been applied to the case of 
semisimple complex Lie algebras in [1]. 

In the next section we implement the above algorithm. In particular, we write down the 
linear system coming from the Jacobi identity and then find the explicit expression for 
c,u and v. We also find the equation, which is the key result of the algorithm, satisfied 
by a. In section 3, as an application, we consider the exponentiation of the Virasoro 
algebra and derive the solution of the corresponding BCH problem. This includes, as 
a particular case, the closed form of the BCH formula for SL 2 (C). We also determine 
three-dimensional subalgebras of the Virasoro algebra satisfying the Van-Brunt and Visser 
condition (11. Ih . In the last section we apply the algorithm to hnd the exponential form 
of SL 2 (C) matrices. Furthermore, we reproduce the same results, using an alternative 
method, and extending them to GL 2 (C) matrices. In this respect, it seems that in the 
literature, similar expressions for 7 G SL 2 (C) are usually given only separately for three 
distinguished cases, depending if + 712721 is negative, vanishing or positive. It turns 
out that the expression of X in 7 = exp(X) has a geometrical representation which is not 
directly evident in 7 . Namely, it turns out that X can be nicely expressed in terms of its 
eigenvalues and of the fixed points z±, solutions of the equation z = jz, where jz is the 
Mobius transformation 


7ii^ + 7 i2 /, , ON 

jz := - . ( 1 - 13 ) 

I21Z + 722 

Such a geometrical representation of 7 may be of interest in the framework of Uniformiza- 
tion theory and Conformal Field Theories. 


2 Implementation of the algorithm 

Let us write down the linear system that follows by the Jacobi identity fll.lip 

uw + mz = 0 , 
vm — wp + n{z — u) = 0 , 
pu -\- zv = 0 , 

c{w + m) + e{z — u) — d{p -|- u) = 0 . (2.1) 

Replacing X and Y on the right hand side of fll.Sp by their expressions in terms of X, Y 
and /, and comparing the result with the direct computation, by fll.51) . fll.jp and fll.lOp . 
of [X, V], yields 

c = {hp{z, w) - gp{z, w)la{u, v)m)c + {ha{u, v) - g^iu, v)lp{z, w)p)d + g^iu, v)gp{z, w)e , 
u = hjs^z, w)u + g/ 3 {z, w)m , 

V = ga{u,v)pYha{u,v)z , 

uK{u, v) + vhf}{z, w) = ga{u, v)hf}{z, w)v + g^iu, v)gf}{z, w)n -F V(u, v)gp{z, w)w . 

(2.2) 
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The first three equations fix c, u and v in terms oi a = 1 — (3. Replacing the expressions 
of u and v in the fourth equation provides the following equation for a 


K{u,v)[hp{z,w){u+z)+gp{z,w){m-w)]+ga{u,v)[hp{z,w){p-v)-gp{z,w)n] = 0 . (2.3) 


This is the basic equation of the algorithm and is further investigated, together with the 
linear system (|2.1|) . in [5]. Note that 


V — au sinh(n/2) 

V sinh[(n — Q;n)/2] 

V — au sinh(Q!-u/2) 

u sinh[(n — aM)/2] 


(2.4) 


We then obtained a closed form of the BCH formula for cases in which the commutator 
contains other elements than the ones in the commutator. The above algorithm can be 
extended to more general cases, e.g. by considering decompositions like the one in (II.dh 
for 

exp(W) • • • exp(X„) . (2.5) 


3 Exponentiating the Virasoro algebra 


In this section we apply the algorithm, leading to closed forms of the BCH formula, in 
the case of the exponentiation of the Virasoro algebra 



[Cj, Ck] = {k- j)£j+fc + ^{k^ - k)6j+k,oI , 

(3.1) 

j, k E Z. 

In particular, we hnd the closed form for W in 



exp(X) exp(V) exp(Z) = exp(lT) , 

(3.2) 

where 

X := X-kC-k , Y := Ao^o , Z := XkCk ■ 

(3.3) 


This is particularly interesting because we do not know alternative ways to get it. The 
case of SL 2 (C) follows straightforwardly since the sl 2 (M) algebra 


[Tj, Z//j] (A; j)Tj_|_fc , 


(3.4) 


j, k = —1, 0,1, is a subalgebra of (13.11) . Note that setting = Li, H = —2Lq and E+ = 
—L_i, reproduces the other standard representation of the sl 2 (IR) algebra [H, E+] = 2E^, 
[E^,E-] = E[ and [H,E-] = —2E-. Note that. 


v2k r 

[X, V] = fcAoX , [V, Z] = k\,Z , [X, Z] = \_k\k + —{e - k) 

LAq iz 


(3.5) 


where, besides c = (i = n = tc = 0, we have, consistently with the Jacobi identity, 
m = p = 0. The other commutator parameters are 


u = z = 
n = X-k 


kXo , 
2k 

Afc — 

Ao 
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e = 


X-kXk—{k^ 


k) . 


(3.6) 
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By fl2.3l) it follows that a satisfies the equation 


nga{u,Q)gp{u,Q) = 2uK{u,Q)hp{u,Q) , 


so that, using 


h„(M,0) = Q!, 5 (o(m,0) = 


au 


1 — e~ 


(3.7) 

(3.8) 


and recalling that A_ := AqO, A+ := Ao/5, one gets 


-kx+ 1 + 6 ~ k^\-k>^k ± a/( 1 + e - PA_fcAfc)2 - 4e 

e — 2 . (3.yj 

Next, observe that, by fl2.2l) . u = fcA+, h = /cA_, and Ck = c = ega{kXo, 0)(7^(/cAo, 0), that 

is 

, (3.10) 


— 


X-X-k A+Afc c 
1 - e-^^- 1 - e-^^+T 2 ' 


that, by fl3.9p . is equivalent to 

X—kXk I A^ 

C-k 


A_ Ac 


A . - A_ V 1 - e-^^+ 1 - e-^^- / 12 


-(r-r) 


(3.11) 


Finally, Eq.([ri]) yields 


exp(A_fc£_fc) exp(Ao£o) exp(AfcA) = 


A+-A 


kX-kC-k + ( 2 — e — e )-^o + kXkCk + Ckl 


(3.12) 


In the case Aq = 0 we have 
exp(A_fe£_A;)exp(AA:£fc) = exp 
with 


A_ 


sinh(/cA^ 


-{kX-kk^-k + k"^ X-kXk^o + kXkCk + Ckl) 


Ck = A-tAt—(A^* - k'^) . 


(3.13) 

(3.14) 


Let us report the relevant case corresponding to SL 2 (C), obtained by setting fc = 1 in the 
above formulas. We have 

exp(A_iL_i)exp(AoLo)exp(AiLi) = 


exp 


A+-A- 

g-A_ _ g-A^ 


X—\L—\ -|-(2 — e — e ^ j+o + AiL 


and 


exp(A_iL_i) exp(AiLi) = exp 


sinh(A^ 


-(A_iL_i + A_iAiLo + AiLi) 


(3.15) 


(3.16) 
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Interestingly, the algorithm for the BCH formula may be extended to any three-dimensional 
subalgebras of the Virasoro algebra. In this respect, it can be easily seen that the highest 
hnite dimensional subalgebras of the Virasoro are the four-dimensional ones generated by 
C-n, Cq, Cn and the central element /, for all n G Z\{0}. Above we solved the BCH 
problem for all such subalgebras. Since all the remanent non-trivial Virasoro subalgebras 
are either three-dimensional, each one containing /, or two-dimensional, it follows that 
all of them satisfy the condition (II.Ih and therefore the corresponding BCH problem of 
hnding Z such that exp(X)exp(V) = exp(Z) is easily solved by (II.2p . In this respect, 
note that it can be easily seen that the two dimensional subalgebras are generated by 
and £o) for all n G Z\{0}. A two-parameter family of three-dimensional subalgebras of 
the Virasoro algebra, is the one where each subalgebra is generated by 

X„(5, 6 ) := 5Co + , X_„(e, 5) , / , (3.17) 

whose commutator is 

[Xn{6, e),X_n(e, 5)] = -neX^iS, e) - n6X_n{e, 6) + - n)I . (3.18) 

Other three-dimensional subalgebras of the Virasoro algebra are the one-parameter family 
of subalgebras, each one generated by 

Xn{o!-) '■= c,2n + OlLn + , V_„(a) := L-n H—TIq , 7 , (3.19) 

9 a 

whose commutator is 

77/2 C 

[Xn{a), Y_n{a)] = -6-Xn{a) - -na‘^Y_n{a) + a—{n^ - n)I . (3.20) 

a 9 12 


4 Geometrical constructions 


Note that, in the case of SL 2 (M), therefore for real A^’s, depending on the values of 
the Afc’s, the factor in 03.151) may take complex values. This is in agreement 

with the non-surjectivity of the exponential map for sl 2 (M) into SL 2 (M). In particular, 
exponentiating sl 2 (M) cannot give SL 2 (M) matrices whose trace is less than —2. The case 
with trace —2 and non-diagonalizable matrices is critical, both for SL 2 (M) and SL 2 (C). 
Let us express Eq. (l3.15p in terms of the associated SL 2 (C) matrix 


7 = 


A B 
C D 


(4.1) 


Replace then the L^’s in the left hand side of (I3.15p by their matrix representation 


L-i = 

using 

exp(L_i) = 


0 -1 

0 0 


Ln — 


0 


, exp(Lo) = 


e 2 


0 

1 

62 


Li = 


0 0 
1 0 


exp(Li) = 


1 0 
1 1 


(4.2) 


. (4.3) 
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Comparing the result with 7 yields 

A = - A_iAi) , B = -A_ie^°/2 , C = , D = , (4.4) 


so that 


e = 


t ih \/B — \ 
D 


(4.5) 


where t := ^try. Since the eigenvalues of 7 , solutions of the characteristic polynomial 
— 2tu + 1 = 0 , are u± = t ± — 1, we have 




(4.6) 


SO that 

exp(A_iL_i) exp(AoIvo) exp(AiLi) = exp 

■ Z/+ — z/_ 

Note that when 7 is parabolic, that is for |f| = 1, and therefore = z/_, we have 

ln(z/+/i/) 


{^^^^^ 7 ^[C'Li + (Z1-AI)Lo-5L_i]} . (4.7) 


U: — V- 


= 2 sgn(f) 


(4.8) 


Also, note that when 7 is elliptic, that is for |f| < 1 , one has = pe*®, so that 

ln(z/+/i/_) 9 


Eq. fi 4 . 7 p implies 


A B 


= exp 


C D 

equivalently, using = 1 , 


Z/+ — p sin 6 

\Yi{v+/v_) ( {A-D)/2 B 


U: — iy_ 


C {D-A)/2 


7 = exp 


ln(f + y/B — 1) 


(7 - th) 


(4.9) 


(4.10) 


(4,11) 


which is indehnitely iterable by replacing 7 on the right hand side, by its exponential 
form, that is by the expression on the right hand side itself. 

The relation fl4.1ip can be derived in an alternative way. First note that 7 = exp(X) 
does not uniquely £x X. For example, exp(X) = exp(X) exp(27rifcJ2) = exp(X + 2'Kikl2), 
fc G Z. However, X can be consistently hxed to be traceless, so that, being proportional 
to I 2 , one has 7 = exp(X) = aX + tl 2 . Therefore 7 = exp[a “^(7 — tl 2 )\ for some a. For 
distinct eigenvalues the diagonalization of both sides (by the same matrix), reproduces 
(14.1111 . since it fixes 


— V- 

a = —;- 

2 In 

A particular case of (I4.1ip is when D = A~^ 


(4.12) 



r / 

1 

2C 

2B 

\ 1 

7 = exp 

In A 

A-A-i 

-1 



L \ 

A-A-^ 

/ J 


(4.13) 
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cannot 


where, since AD — BC = 1, either i? = 0 or C = 0. It follows, that when D = A~^, 

with = 1, 7 admits exponentiation only if A = 1. For example, ^ ^ 

be expressed as the exponential neither of sl 2 (M) nor of sl 2 (C). Of course, this is not a 
problem in the case such a matrix is seen as an element of PSL 2 (M) = SL 2 (M)/{±I}. 

Let us derive a more geometrical representation of (Id.llh . useful, e.g., in the framework 
of Uniformization Theory and in Conformal Field Theories. Consider the Mobius trans¬ 
formation 

Az + B 

: = 


Cz + D' 

and the solutions of the hxed point equation '-)z = z 


z± = 


[A - D)l2±y/W^ 

C 


If 5 = 0, then 


f dl 

0 \ 


nn(i/+/n_) 

f z+ + Z- 0 \ 


0) 

= exp 

z+ — 

V 1 ). 


(4.14) 


(4.15) 


(4.16) 


otherwise, in the case of Mobius transformations, one can fix 5 = 1 to get 

■ln(z/+/z/_) 


A 1 
C D 


= exp 




2 :+ + 
1 


2 |.— 2 : 


— 


(4.17) 


We conclude by observing that our hndings trivially extend to GL 2 (C). In particular, 
multiplying both sides of fld.lip by IfI := defy, one gets for 7 G GL 2 (C) 


7 = exp 


ln(f ^/W 




—(7 - 1/2) + bn(h -|)/2 


(4.18) 


Acknowledgements 

It is a pleasure to thank Pieralberto Marchetti, Leonardo Pagani, Paolo Pasti, Dmitri 
Sorokin and Roberto Volpato for interesting discussions. 


8 






















References 

[1] A. Van-Brunt and M. Visser, arXiv: 1501.02506 

[2] A. Van-Brunt and M. Visser, arXiv: 1501.05034 

[3] M. Matone, arXiv: 1503.08198, 

[4] M. Matone, arXiv:1504.05174, 


9 


